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Abstract. We investigate the variations of broadband luminosities of low and intermediate 
, redshift galaxies due to the stochastic nature of the neutral hydrogen distribution present in 

form of Lya absorbers in the intergalactic medium. This effect is caused by variations in the 
distribution and properties (redshift, column density, Doppler parameter) of the absorbers along 
different lines-of-sight out to a given redshift. Using a set of observationally constrained redshift-, 
column density- and Doppler parameter distributions we perform Monte Carlo simulations for a 
large number of lines-of-sight towards galaxies at a given redshift z em and calculate attenuated 
FUV/NUV magnitudes and corresponding 1-, 2- and 3<j variations. We predict significant 
variations in luminosity ranging from -0.22 to +0.48 magnitudes at the la level for 
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1^-^ , galaxies at z em — 1.5 in GALEX data. 
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■ 1. Introduction 

The observed photometric properties of local and distant galaxies are important to 
constrain evolutionary synthesis models describing galaxy formation and evolution. But 
+-> ■ estimating colors and luminosities of galaxies from models and comparing them to the 

observed ones is a difficult task since, among other effects, the light of a galaxy is partially 
absorbed by the intervening neutral hydrogen (Hi) present in the intergalactic medium 
(IGM) in form of discrete Lya absorbers. The distribution of these Lya absorbers along 
the line-of-sight (LOS) is stochastic in nature, i.e. the population of Lya absorbers char- 
acterized by their redshift £ a bs> column density Nm and Doppler parameter 6 is unique for 
each LOS. Since the observation of individual galaxies implies necessarily the observation 
along different lines-of-sight, we expect that the stochastic nature of the Lya absorbers 
distribution, especially of those at the highest column densities, should cause a significant 
scatter in the observed broadband colors, even for galaxies with similar intrinsic colors. 
Furthermore, we expect this stochastic effects to be strongest along the shortest LOSs 
and, hence, to be only observable in the UV as e.g. with GALEX. 

In order to correctly account for the absorption of light caused by this intervening ab- 
sorbing Hi-systems, a detailed knowledge of their properties (i.e. z a i, s , Nm, b) is required. 
It has been found from statistical analyses of quasar absorption lines that, when averaged 
over a large number of LOSs, the mean distribution and properties of Lya absorbers are 
well described by a redshift- and column density-dependent power law of the form of eq. 



12.11 (see e.g. Kim et al. (1997) and references therein). Using different sets of parameters 
for this power law derived from observations, there have been various studies (see e.g. 
M0ller fc Jakobsen (1990)1 |Madau (1995)1 |Bershady, Charlton fc Geoffrey (1999)1 here 



after BCG99), which, with a different purposes and approaches, have made an effort to 
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account for the mean attenuation of the light of distant galaxies due to intervening Lya 
absorbing systems and the dispersion around this mean. In this work, we make a further 
effort to calculate the stochastic effect of the Lya absorbers on the photometric proper- 



ties of galaxies. Applying some of the m ethods develope d by M0ller fe Jakobsen (1990) 
and BCG99, we recover the formalism of Madau (1995)[ i.e. we model the mean attenu- 
ation and any desired a level around this mean as a function of redshift. The advantage 
of our approach is that it is of general use in that it is completely independent of the 
galaxy input spectrum, and that it includes the corrections concerning the calculation of 
the a scatter pointed out by BCG99. As an application of our model, we estimate mean 
magnitudes and scatter at the 1-, 2- and 3<r level in the GALEX FUV/NUV broadband 
filters for galaxies with flat input spectra at 0.3 < z em < 1.5. 



2. The Model 

We present Monte Carlo simulations for the distribution of intervening Hi-systems 
along random LOSs out to a given redshift z em (see also M0ller fc Jakobsen (1990)[ 



BCG99). As input we use differential distributions for the number density of clouds as a 
function of redshift and their column densities of the form 

and a redshift-dependent Gaussian distribution for the Doppler parameters, all derived 



from observations (Kim et al. (1997) BCG99). We generate an ensemble of 4 • 10 3 LOSs 
out to a given z em , each of them containing a random absorber population drawn from the 
input distributions. For each absorber along each LOS and we calculate the photoelectric 
and Lyman-Series absorption assuming Voigt- profiles and a flat input spectrum with rest- 
frame resolution AA = 0.5 A. We obtain a set of transmission values at each point of our 
wavelength baseline and get in this way an ensemble of transmission values in the range 
[0, 1] as a function of wavelength. 

On this basis we define the Transmission Function T(e _T | A ; z em ) as the residual 
transmission along a random LOS out to z em at a given wavelength A as a function of 
the absorption coefficient e _r . It gives, up to a normalization constant, the probability 
for the absorption coefficient e~ r to take a certain value in the range [0, 1] at a given 
wavelength along a random LOS out to z em . Since there is no known analytic expression 
for T(eT T | A ; z em ), we determine it numerically from our simulations. It turns out that 
this distribution function is rather peculiar, and the calculation of relevant statistical 
moments for the absorption assuming Gaussian statistics may lead to incorrect estimates, 
as already noted by BCG99. 

Let us define Tq as the q transmission quantile via the integral equation 

q = C Q [ q T(l-r | A;z em )dr (2.2) 



where Co is a normalization constant, and the mean value of T(e T | A ; z em ) as T50; and 
the 1-, 2- and 3er ranges around this value as T5o±34.i3, T5o±47.72, and t5q±49.86; respec- 



tively. Using this definition we calculate mean attenuation curves (see Madau (1995) 1 and 
the corresponding 1-, 2-, and 3cr levels for any desired redshift, which can be applied to a 
given spectrum to estimate the amount of absorption due to the stochastic distribution 
of the Lya absorbers and their properties (see figQ . 
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Figure 1. Attenuation curves for z em = 1.5: mean and 1-, 2- and 3a variations (see ea. 12.21 . 
For clarity, the curve corresponding to 'mean - 3<r' is not shown. 



z em Passband mean la 2a 3a 

0.3 FUV -18.82 

0.3 NUV -20.26 

0.5 FUV -18.69 

0.5 NUV -20.25 

1.0 NUV -20.25 

1.5 NUV -20.02 



Table 1. Mean and 1-, 2- and 3a variations in luminosity (Mab) in the GALEX FUV and NUV 
passbands for a flat input spectrum at four different redshifts z em £ {0.3, 0.5, 1.0, 1.5}. In the 
cases where no value is given, either there is no significant variation or the flux inside the given 
filter has dropped to zero and the corresponding luminosity is not longer defined, being the case 
when the Lyman-break has already cross the filter or due to strong line absorption in the given 
filter. 

3. Preliminary Results & Conclusions 

We apply our method to calculate the effect of the Lya absorbers on the FUV/NUV- 
broadband luminosities of a galaxy at redshifts z e {0.3,0.5,1.0,1.5}. Our results are 
given in Table 2. A careful inspection of these results shows that the contingent presence 
of Lya clouds with the highest column densities can cause substantial variations in 
luminosity already at the la level, ranging from -0.22 to +0.48 magnitudes. 
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